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ABSTRACT 

I consider how physical processes scale over eight orders of magnitude in black 
hole mass, from stellar masses in gamma-ray bursts (GRB) and black-hole X-ray 
binaries (BHXRB) to supermassive active galactic nuclei (AGN). Accretion rates 
onto stellar mass black holes range over more than sixteen orders of magnitude, 
from the lower luminosity BHXRB to GRB. These enormous parameter ranges 
correspond to qualitative as well as quantitative differences in behavior. The 
fundamental questions involve the balance between nonequilibrium and thermal- 
ized plasmas. When energy fluxes exceed a critical value ~ 10^^ erg/cm^s, as in 
GRB, a black-body equilibrium pair plasma forms. At the lower fluxes found in 
AGN, BHXRB and microquasars, accretion power electrodynamically accelerates 
a small number of very energetic particles, explaining their non-thermal spectra 
and the high energy gamma-ray emission of blazars. Ultra-high energy cosmic 
rays may be accelerated by massive black holes, otherwise undetectable, with 
very low thermal luminosities. New-born fast high-field pulsars may be in the 
black-body equilibrium regime, resembling SGR in permanent outburst. I also 
consider the question, significant for the acceleration of nonthermal particles in 
GRB outflows, of whether collisionless plasmas interpenetrate rather than form- 
ing hydrodynamic shocks, and propose this as an alternative to internal shock 
models of GRB. 

Subject headings: Accretion: Accretion Disks — Galaxies: Active — Galaxies: 
Nuclei — Gamma Rays: Theory — Gamma Rays: Bursts — Acceleration of 
Particles 
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1. Introduction 

Hundreds of models of active galactic nuclei (AGN) and quasars have been published 



since the pioneering work of (ISalpeter 19641 ). yet none of them is completely satisfactory and 
generally accepted. This papeil^investigates the general scaling laws that connect the very 
different parameter regimes appropriate to AGN and to BHXRB, including microquasars. 

GRB, AGN and BHXRB have certain qualitative similarities, despite enormous quan- 
titative differences in their luminosities, masses, and time scales. They all are powered by 
accretion onto a central black hole. They all produce nonthermal radiation. GRB and AGN 
and some BHXRB show evidence for relativistic bulk motion. They all fluctuate irregularly 
in intensity. AGN and some BHXRB (microquasars) are directly observed to produce narrow 
jets, while energetic arguments and modeling of their afterglows by (d) indicate that GRB 
are also strongly coUimated. 

These similarities suggest that it may be useful to investigate how a common accretion 
disc model may manifest itself differently in these objects. The present discussion builds 
upon the (0; |29| ) model for AGN, although these authors assumed (probably unnecessarily, 
except to make the problem tractable) that the magnetic dipole moment is aligned with a 
rotational axis. These models may be scaled to stellar mass BHXRB such as Cyg X-1 and 
the microquasars. 

A single parameter divides all accretion discs into two classes. The fundamental differ- 
ence between these classes of models is that in GRB and in SGR in outburst (and probably 
also in young fast high-field pulsars, as yet unobserved) the relativistic wind is thermal- 
ized to an equilibrium pair plasma, while at the lower power densities of AGN, BHXRB 
and observed radio pulsars the wind remains transparent and very far from thermodynamic 
equilibrium with the radiation field, as a comparatively few particles are accelerated but to 



high energies. The condition for thermalization (jl8l : l2d : l22l : l23l ) is that 

I = I/Ichar > 1, (1) 

where I is the energy flux (power emitted per unit area, whether carried by particles, ener- 
getic photons, or low frequency or even DC Poynting flux), so that / exceeds a characteristic 
value 



'char = CTSB 



nipC^ 1 4 X lO^^erg/cm^s 

^ (1 + 0.05 In M/M©)^' 



_A;Bln(6(M/M0)iVB0«G«2 

where asB is the Stefan-Boltzmann constant, Nbq ~ 1.2 x 10^'^ is the number of baryons 
in a Solar mass, ac = Gml/hc ~ 1.76 x 10~^^ is the "gravitational fine structure constant" 



^The present paper is the "forthcoming work" cited as Katz (1997b) by (j23r ). 
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defined for tlie electron mass. Tliis corresponds to an equivalent black body temperature 
exceeding tlie cliaracteristic value 

2.9 X 10^°K 
1 + 0.05 In M/Mq 

at which a black body equilibrium pair plasma is opaque to Thomson scattering over the 
characteristic length scale (three Schwarzschild radii) of the source. Because of the steep 
Boltzmann factor in the equilibrium pair density, the characteristic temperature and intensity 
are almost universal constants, only weakly dependent on the mass of the source, even over 
the range from IMq to IQi^Mq. 

The thermalization of energetic particles by scattering in regions of high energy density is 
auto-catalytic: scattering increases the number density of potential scatterers (by processes 
such as double Compton scattering, pair production, and curvature radiation). However, 
these processes only run away to full thermalization if X ^ 1. 

In these models of AGN and BHXRB the electromagnetic energy is converted to the 
energy of accelerated particles close to the black hole. This is in contrast to GRB in which 



electromagnetic energy is converted to a black-body equilibrium pair plasma (1351 ) 



In another class of models of AGN and BHXRB the disc radiates vacuum electromag- 
netic waves instead of energetic particles. At much greater radii these waves accelerate 
particles, just as in GRB the pair plasma accelerates particles in distant collective interac- 
tions. These two classes of models can be comparably efficient particle accelerators. I do 
not consider the vacuum wave model further because it is less closely analogous to the GRB 
model (which cannot be a vacuum wave model because the energy density leads to creation 
of an equilibrium pair plasma), and because external plasma injection or pair breakdown are 
likely to fill the wave zone with energetic particles. 

The remainder of this paper is chiefiy concerned with scaling of a variety of physical 
processes over the very wide range of parameters encountered. It would be difficult to 
say anything qualitatively new about any of the models of these objects, which have been 
extensively developed over 45 years, but it may be useful to re-examine their scaling. To 
this end it is necessary to review the fundamental physical processes involved. 

There are also temporal similarities among the various classes of black hole accretion 



discs. In one BHXRB (Cyg X-1; (1431)) a non-zero time skewness was measured in an X-ray 



time series. Time skewness of the same sense is found ((1321)) in many GRB. Searches for 
time skewness in AGN time series have so far been unsuccessful ((j36l)). but do not exclude it. 



In fact, it was once suggested ( (1331 : |3J)) that Cyg X-1 might be a source of GRB. Although 



such Galactic stellar mass X-ray sources are now known not to be the origin of GRB, Cyg 
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X-1 does show outbursts and flaring behavior in comparable energy bands {(^^] 



2. Homopolar Generators 

All of the objects we discuss contain rotating, electrically conducting, magnetized mat- 
ter, and therefore are homopolar generators with potentials Vq ~ J{v/c)Bdr. In general, 
such generators produce a power P = Vq /Z, where the load impedance Z = E / H = Zj. + iZi 
is written as the sum of reactive {Zi) and dissipative {Z^) parts. A fundamental assumption 
of nearly all astrophysical models of such objects is that \Z\ ^ Zq = Att/c = 4.19 x 10^^'^ 
sec/cm, the impedance of free space in c.g.s. units {Zq = ^/Jlo/eo = 377 in m.k.s. units). 

The detailed or microscopic justification of \Z\ ~ Zq is model-dependent. A similar 
result will always be found if the particle energies are large enough that particle multipli- 
cation (by processes such as double Compton scattering and collisional and radiative pair 
production) turn the insulating vacuum into a conducting relativistic plasma. Zi ~ Zq is a 
general result obtainable from Maxwell's equations by estimating |V| ~ d/d{ct). As long 
as the system size exceeds cAt, where At is a characteristic (rotation or fluctuation) time 
scale, the reactive impedance Zi ~ Zq will always be in series with Zr. 

The dissipative impedance Z^ may be small in a relativistic plasma, but we hypothesize 
that when the fields are large particle acceleration and multiplication will generally lead to 
Zr ~ Zq. This hypothesis, essentially dimensional, is central to the model of this paper and 
is implicit in most pulsar and AGN models 

(SiEi)). 

For example, application to pulsars 

gives 

P r-u {Ery [ -T B ] ~ ~ (4) 

Z An An \c / Anc Ancr^ Anc^ 

where /z is the magnetic moment and S ~ /i/r^ is evaluated at the radius of light cylinder 
r = c/uj. This is comparable both to the expression for magnetic dipole radiation of a 



misaligned rotor in vacuum and to the result of (jl3l ) for an aligned rotor when the vacuum 
fills with a relativistic plasma. It is also unavoidable on dimensional grounds, provided that 
the fields are large enough to drive these processes. 

In the laboratory Z^. 3> Zq when the vacuum or materials do not break down into a 
conducting plasma, as when fields and potentials are small and insulators remain insulators, 
and frequently ^ Zq when conductors are not relativistic. In relativistic astrophysics 
Zr ~ Zq follows from assuming that the power carried to r > cAt by the Poynting vector 
is dissipated somewhere, rather than remaining forever as electromagnetic field energy or 
returning (like an unphysical advanced time solution) to the source. 
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The first nontrivial part of tliis problem is to estimate B, wliicli can be done on tlie 
basis of dimensional arguments from tlie total disc power. The second, and most critical 
(but generally unsolved), part is to determine how the potential drop is distributed in space 
and how the dissipated power is apportioned among the accelerated particles. 



3. Nonthermal efficiency 

All three classes (AGN, BHXRB and GRB) of objects show a significant amount of 
nonthermal emission. The nonthermal efficiency, which may be defined as the fraction of the 
emitted power appearing as non-Planckian radiation or as particle distributions that are far 
from thermal equilibrium (typically power laws over orders of magnitude in energy) is 0{1). 
In GRB the observed emission appears to be nonthermal, although it is not the primary 
radiation emitted by the central engine but rather the consequence of particles accelerated 
in the relativistic outflow; the high energy density and intensity at the source thermalizes 
the relativistic wind. Still, the wind is generally believed to be produced by a fundamentally 
nonthermal process, the coherent radiation of electromagnetic energy by the fields of the 
central engine. In many AGN a substantial fraction of the power appears as nonthermal 
visible synchrotron radiation or high energy gamma-rays (I28!). 

The case for nonthermal X-ray emission in BHXRB is plausible but less compelling. 
It is a natural explanation of the complex multi-component spectra frequently observed, 
but it may also be possible to explain such spectra as the sum of thermal spectra emitted 



by matter distributed over a range of temperatures. For example, ( l30l ) gives the classic 
example of an apparently non-thermal power law spectrum that results from thermal emission 
from an accretion disc whose temperature is a power law function of radius. Only if the 
emission extends to very high energy, or is associated with low frequency emission of very 
high brightness temperature (as in radio pulsars), can the case for nonthermal processes be 
considered compelling. 

As examples of such a compelling case, the superluminal radio components and jets 
present in some BHXRB (microquasars) certainly require acceleration of relativistic particles. 
Some, such as Cyg X-3, also show strong outbursts of nonthermal radio emission (there is no 
direct evidence this particular object contains a black hole, and qualitatively similar effects 
may be produced by accretion onto a low magnetic field neutron star), but this radio emission 
is only a very small fraction of their power output. 

In GRB the electrodynamic efficiency eg (Equation 5 of (23)) is not directly measured 
because virtually all the thermal radiation emerges as neutrinos (the thermal luminosity 
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is subject to an Eddington limit that is ~ lO^^M/M© erg/s for electromagnetic radiation 
but many orders of magnitude greater for neutrinos because of the tiny neutrino opacity 
of matter), and is essentially undetectable. For the lower energy density accretion flows of 
AGN and BHXRB neutrino emission is negligible, and the thermal radiation produced by 
viscous heating is directly observable. 

(23) argued that eg ~ 0.1-0.5 is likely, independent of the magnitude of the magnetic 
field (but depending on its unknown orientation and spatial structure). This is consistent 
with observations of AGN; the likelihood of relativistic beaming precludes quantitative com- 
parisons. This range of eg is also consistent with observations of Cyg X-1 and other BHXRB 
if the harder components of their spectra are either nonthermal or the thermal emission of 
optically thin matter heated by nonthermal particles. 

The measured nonthermal efficiency is also affected by radiation trapping ((ll7i)). As is 
well known (and the subject of an extensive literature) if the mass accretion rate exceeds 
the nominal Eddington rate the excess mass is readily swallowed by the black hole, but 
the luminosity Lth emergent in thermal radiation that diffuses through the accretion flow is 
limited to slightly less than the Eddington limiting luminosity Le {cf. (5)). An analogous 
limit involving the neutrino Eddington limit applies to the unobserved neutrino luminosity. 
In AGN and BHXRB there can be an apparent nonthermal efficiency eg — > 1 and Prw ^ Lth, 
because the nonthermal wind luminosity Prw is proportional to the accretion rate and is not 
subject to the Eddington limit, while the emergent thermal luminosity Lfh cannot exceed 
Le even if the viscous power dissipated Pvisc ^ Le- 



4. Particle acceleration 

In this section I collect and summarize results (many of which are spread throughout 
the very extensive literature) concerning nonthermal acceleration processes that may occur 
in black hole accretion discs. The purpose is to develop criteria to determine if the condition 
1^1 necessarily leads to black-body equilibration, which is energetically allowed, but not 
required by any fundamental physical law. Unfortunately, the rates of many processes depend 
sensitively on the angular, spatial and (in the Klein-Nishina case) spectral distributions of 
any soft photon background flux, such as may be emitted by clouds of comparatively cool 
thermal matter. These fluxes depend on detailed modeling that cannot be developed from 
flrst principles, and for which we do not generally have sufficient information to enable even 
phenomenological models. 

Pulsars and GRB (in the present model) have large electric fields which lead to pair 
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production. A supermassive black hole in a galactic nucleus, or a stellar mass black hole in 
a mass transfer binary, is surrounded by a complex accretional gas flow. Significant sources 
of mass include the companion star (in the BHXRB), the galactic interstellar medium and 
possibly disrupted stars (in the AGN) and the surface of the outer parts of the accretion 
disc. Although the fiow is not understood in detail, it is plausible that some of this gas has 
sufficiently little angular momentum (or loses its angular momentum at large enough radii) 
to permit accretion on the axis of rotation, and can fill all directions around the black hole, 
as was found in the calculations of (5). A pulsar-like vacuum is not likely. The space charge 



density required to neutralize the corotational electric field ((]l3|)) is small, and may readily 
be supplied by this plasma. Pair production is therefore not required for the extraction of 
energy from the rotating disc in AGN and BHXRB, although it may occur. 

The homopolar generator of a rotating magnetized fluid ( (fisl : 0; [29I )) implies a large 
electric field with a nonzero curl in an inertial (observer's) frame. This cannot be canceled 
or screened everywhere by space charge density. In some regions the electric field must be 
nonzero (and large); these regions are conventionally called "gaps" because the field purges 



them of plasma. Charged particles entering or created ( (l4ll : l38l) ) in these gaps are accelerated. 
The fundamental unsolved problems of compact magnetized rotating astronomical objects 
are the distribution of the homopolar potential drop and the acceleration of particles in gaps. 

If the entire circuit were composed of dense thermal plasma then Zr would be the plasma 
resistivity, which would be -C Zq. Zi would depend on the circuit dimensions. If these were 
much less than the characteristic light travel dimension for time variations (as often the 
case in the laboratory, where circuits of cm or m dimensions may oscillate at 60 Hz, or 
even be nearly DG) then Zi may be tiny (^ IQ). The power P = V"g^Z,./|Zp ^ VQ/Zr 
could be extremely high. This may describe some GRB (perhaps short GRB or subpulses 
in longer GRB, though it is hard to be sure because their magnetic fields are unknown), 
in which vacuum gaps may fill with equilibrium pair plasma, or electromagnetically driven 
supernovae. However, if the circuits close at the light travel distance or are open beyond 



it, as in classical pulsar models (Il3l ). then Zi ~ Zq and P < Vq/Zq (as in Eq. H]). This is 
probably the case for AGN and BHXRS, in which Vq is largely determined by fields near the 
light travel distance. 

Near a luminous object accelerated electrons and positrons are slowed by Gompton 



scattering on the thermal radiation field ( (jl6l : l25l)). As a result, they are not accelerated 



to the limiting energy Eq = cVq ~ e J {v/c)B dr. As an extreme bound we assume that 



this electric field is parallel to B. In the presence of a significant ((jl3l)) density of charged 
particles the parallel component of E will be much less. The following results will remain as 
upper bounds and may be useful as such, even though they are likely to be overestimates of 
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the actual energies achieved. 



4.1. Electron Retardation by Compton Scattering; Thomson Limit 

If the Thomson cross-section is applicable the energy loss length of an electron with 
Lorentz factor 7 in an isotropic radiation field of intensity Lth/ {^Tir'^) near a mass M is 

rriel Le rc^ 
nip 7 Lth GM 

The first two factors are each <^ 1, and the next two are not much greater than unity in the 
inner disc of a luminous object, so that fc" <C r. This result is also approximately valid for 
anisotropic radiation fields, except in the extreme case of a particle moving accurately in the 
direction of a narrowly collimated beam of radiation. 

Given a value for the magnetic field B, it is possible to calculate the maximum energy 
an electron achieves, and to estimate its radiation. To estimate B equate the magnetic stress 
to that required to supply a bolometric luminosity Lb (including relativistic wind, thermal 
radiation and radiation advected into the black hole) yields 

52 1 Lt, /GMy/' r , U / lOGM y/' r 10«Mq erg 

87T^2LE\rc^) hGMn^ ^ Le \ rc^ J lOh M cm^' ^' 

where k is the Thomson scattering opacity and h the disc thickness. We assume only that 
the viscosity is magnetic and that B^ ~ {Bj.B^). The magnetic field does not depend on 
an assumption of equipartition or on the value of a, and is derived from the accretion rate 
implied by Li, alone. This estimate of B is supported by polarimetric observation of the 



BXHS Cyg X-1 (djj)) 



Although the estimate of B (and especially of (BrB^)) is robust, the process by which 
accretional energy is dissipated remains mysterious. It is plausible, or even likely (in analogy 
with the magnetic heating of the Solar convective zone or the tidal dissipation in the ter- 
restrial oceans), that dissipation occurs in a disc corona rather than in its deep interior. In 
a low density corona the power dissipated may be coupled to a few particles accelerated to 
high energy while power dissipated in the deep interior heats matter in thermal equilibrium. 
Depending on the (unknown) spatial distribution of impedance and magnetic field, these 
loads combine in a complex and unpredictable manner. For example, dense thermal plasma 
of low impedance may be in series with or in parallel to vacuum gaps; in the former case 
most of the dissipation occurs in the gaps, while in the latter case most of it occurs in the 
thermal plasma. 
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The tendency of large electric fields to purge space of plasma, either by collisional resis- 
tive heating and subsequent expansion or by acceleration of coUisionless particles, combined 
with the pinch effect that concentrates the current into narrow filaments, suggests that there 
will be regions of near- vacuum in which a few particles are accelerated to high energy. Hence 
we must consider the limits on particle acceleration in those regions. 

Equating the energy eEic gained in a length ic to the energy 'ynieC^ lost by Compton 
scattering, using Eq. O yields the maximum electron Lorentz factor 



eEr rc^ Li 



1/2 



Using E ^vB/c^ B{GM/rc^ fl'^ and Eq. [6] yields 



e2 ) \Le 



(lo^) (m) (zi. 



where Nb ~ 1.2 x 10^'^ M/Mq is the ratio of the black hole mass to the proton mass. It is 
possible to express 7c in terms of fundamental constants, dropping factors of order unity 
which depend on the properties of the individual object, and noting that if M equals the 
Chandrasekhar mass Men-, then Nb ~ {he/ Gm'^Y^'^: 



4- 1.1. Compton gamma-rays 

Eq. [8] directly gives an upper bound on the Compton scattered photon energy 7cmeC^. 
This limit is ~ 10^^ eV for typical AGN masses and ~ 10^° eV for typical BHXRB masses if 
the factors in parentheses other than M/Mq are of order unity. In AGN typically Lf^ <^ Le 
so that significantly more energetic electrons and Compton scattered gamma-rays may be 
produced. This explains, at least qualitatively, the production of TeV gamma-rays in AGN 
such as Mrk 421 (0)) and Mrk 501 ((37)). 



The actual spectral cutoff depends on the spectrum of thermal photons; if their energy 
huth ~ meC^/jc the cutoff will be ~ 7cw^eC^; if hwth < nieC^/'yc the cutoff will be ~ •y'^huth] 
if fkJth > meC^ lie the cutoff is ~ •yc^eC^ but Eq. [8] then underestimates 7c because of the 
reduction in Klein-Nishina cross-section and the discreteness of the energy loss. The observed 
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spectra of AGN and BHXRB are so complicated (it is also unclear which components are 
emitted at the small radii at which electron acceleration is assumed to occur) that it is 
difficult to be quantitative. Visible radiation from AGN and soft X-rays from BHXRB 
place Compton scattering marginally in the Klein-Nishina range (note that both the black 
body Jiiuth and nieC^ /'jc scale cx M~^^'^), but quantitative estimates depend on the uncertain 
factors in parentheses in Eq [HI 

Because the Compton scattering power of an electron and the frequency of the scattered 
photon each are proportional to 7^, while the rate of energy gain in an electric field is 
independent of 7, there is expected to be a broad peak in uF^, around the cutoff frequency, 
with F^, oc z/^/^ at lower frequencies (the same slope as for synchrotron radiation discussed in 
§4.1.2, for the same reasons). In principle, the weak dependence of 7^ on M, Lfh and Lh in 
Eq. [HI could be tested by comparing AGN and BHXRB if M could be estimated independently 
of the luminosity. 

The majority of the power which goes into lepton acceleration may appear as Compton 
scattered gamma-rays of energy ~ 7(7meC^. This is half Pm, in a proton-electron wind and 
all of Prw if pairs are accelerated. As discussed in §3, this can far exceed Lth if the disc is 
undergoing highly supercritical accretion. This may explain the dominance of the emitted 
power by energetic gamma-rays in some AGN. Supercritical accretion by black holes of 
comparatively low mass also permits more rapid variability than accretion at the Eddington 
limit by a more massive black hole, so that, in principle, masses could be estimated from the 
variability time scale. 



4- 1-2. Synchrotron radiation 

Electrons with Lorentz factors up to that given by Eq. [8] radiate synchrotron radiation in 
the magnetic field Eq. [61 Assuming an isotropic distribution of pitch angles, the characteristic 
synchrotron frequencies extend up to 



. The usual relation between the particle distribution function and the synchrotron spectral 
index for v < Vsynch would predict a spectral index of —1/2, because a uniform accelerating 
electric field produces a particle distribution function with energy exponent zero. However, 
the low frequency component of the synchrotron radiation function increases the spectral 



index to —1/3, in analogy to the predicted ((l2ll)) and observed ( ([ill : [8|; [l|)) low frequency 
spectra of gamma-ray bursts. 
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The ratio of synchrotron to Compton scattering powers, assuming an isotropic electron 
distribution, is the ratio of the magnetic to the thermal energy densities UB/Uth'- 

Psynch ^EB^lL^r^ f GM\ 

Pcompt Uth 2Lthh \ rc^ J 

This is generally > 1. 

Unlike the Compton cutoff, the synchrotron cutoff Eq. [10] is independent of the mass of 
the black hole, and is in that sense "universal", although it depends on other parameters, 
most notably the ratio Lh/Lth- Eqs. [TUl and [TT] suggest the possibility of synchrotron radi- 
ation with significant power up to ~ GeV energies. This is probably a great overestimate, 
both of the power and of the radiation frequency, because the assumption of an isotropic 
distribution of electron momenta is unlikely to be valid. Electrons may be effectively accel- 
erated only parallel to the magnetic field — a component of E perpendicular to B does not 
effectively accelerate charged particles unless it varies at their cyclotron frequency, unlike 
the nearly steady corotational electric field — and readily lose their transverse momentum by 
synchrotron radiation. 

Internal dissipation in the relativistic wind may be as essential to radiation in AGN as in 
GRB, for in them plasma turbulence may partially isotropize the electron distribution, mak- 
ing effective synchrotron radiation possible. If the pitch angles remain small the frequency of 
the synchrotron radiation is reduced and it is emitted nearly parallel to the direction of the 
electrons' motion, the magnetic field, and the Compton scattered gamma-rays. This can be 
described as relativistic bulk motion of the electrons and their associated radiation field, and 
may be necessary to avoid absorption of the photons by gamma-gamma pair production. 



4- 1.3. Curvature radiation 

The electrons also radiate curvature radiation (on the magnetic field lines of radii of 
curvature ~ r) at frequencies up to 

1_ (GM\"' (Uh\-"' (uV" fry/' ruTc^ 
'''''''' ^ 2Arc') \Le) \Le) \h) m'J'ey^iGMy/^ 

Curvature radiation is insignificant, and its power is small, when the electron energy is 
limited by Compton scattering in the Thomson limit. It cannot cause pair production, in 
contrast to the case of pulsars (which have much larger fields than do accretion discs). 
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4-1 .4- Pair production 

The most energetic gamma-rays (produced when the electrons Compton scatter thermal 
photons) of energy ~ 7cmeC^ may produce electron-positron pairs by interacting with the 
thermal photons. The condition for this to occur (assuming an isotropic thermal radiation 
field) 

E^huth ~ 'ycrrieC^huth > {meC^f , (13) 

is equivalent to the condition for the breakdown of the Thomson approximation to Compton 
scattering. It appears to be met in AGN and BHXRB, taking the observed thermal spectra 
and assuming all factors in parentheses in Eq. [HI except that involving the mass, are 0(1). 
If the thermal radiation field is assumed to be a black body then Eq. [T3] can be rewritten 

It is clear that this condition is generally met; if (as is likely) the thermal spectrum is harder 
than that of a black body at the effective temperature the inequality holds even more strongly. 
Pair production by interaction between gamma-rays produced by Compton scattering of the 
accelerated electrons and thermal photons takes the place of pair production by curvature 
radiation which occurs in pulsars. 



4.2. Electron Retardation by Compton Scattering; Klein-Nishina Case 

As pointed out in §4.1.1, Compton scattering of the accelerated electrons by thermal 
photons may be at sufficiently high energy that the Thomson cross-section is inapplicable, 
and the full Klein-Nishina cross-section must be used instead. In this case ic depends on the 
frequency distribution of the thermal radiation as well as on its luminosity. Eq. [S] is replaced 
by 

^ 7 hPth Le huthT ^ 

3\n {-fhuth/rrieC^) rrieC^ LthGMnip 

The dependence on the thermal photon frequency is nearly quadratic, one power coming 
from the reciprocal relation between the photon number density and their frequency (at fixed 
Lth) and the other power coming from the energy dependence of the Klein-Nishina cross- 
section. We have assumed a thermal radiation intensity Lth/4:7ir^, as would be produced 
by energy dissipated at radii comparable to that of particle acceleration. Thermal radiation 
produced at greater radii (for example, in the broad line regions of AGN) is diluted, giving 
a correspondingly lower effective value of Lt^- 
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Because of this sensitivity to uth, information about the spectrum is required to evaluate 
Eq. [151 For example, for near-Eddington limited accretion onto a 10^ Mq black hole ixN ^ ^ 
if huth = 10 KeV, but Ekn r if hvth = 100 eV. A small amount of soft radiation has a large 
effect. For electrons of Lorentz factor ~ 10^, suggested in §4.1, even scattering by visible 
light (below the black body spectral peak for luminous accretion onto a supermassive black 
hole) or the "blue bump" observed in some AGN spectra is in the Klein- Nishina regime, and 
the efficacy of Compton scattering in limiting electron acceleration is sensitive to the actual 
spectral distribution of the radiation. 

The rate of electron energy loss is essentially proportional to the scattering rate, because 
an electron loses most of its energy in a single scattering. The Compton scattering rate is 
given by J aKN{.Fy/v) dv, where axN oc Inijiv /mf,c^) /v. If the thermal spectrum is self- 
absorbed below a frequency Vabsi follows a thin bremsstrahlung or synchrotron spectrum at 
higher frequencies, and Vahs > J^kn where the characteristic Klein-Nishina frequency ukn = 
mgC^/'yh, then most of the Compton energy loss is attributable to photons of frequency 
^ ~ ^abs and Lth in Eqs. |5l [7] and [8] should be replaced by the luminosity at frequencies 
of this order. If Uabs < ^kn then (provided the spectral index in the optically thin regime 
does not exceed 1) most of the Compton energy loss is attributable to photons of frequency 
v ~ vrn^ and Lt^ should be replaced by the luminosity at frequencies of this order. These 
are reasonable rough approximations to the full integrals of the Compton energy transfer 
function ((i3f)) over the radiation spectrum. 



4-2.1. Ultimate (Goldreich- Julian) limit 

The reduction in the Compton energy losses when the bulk of the thermal luminosity 
is emitted at frequencies above ukn implies that the maximum Lorentz factor is larger than 



indicated in Eqs. O [8] and [91 The available potential drop is given by (Il3l ) and accelerated 
particles of mass m reach a limiting Lorentz factor 

iGJ-^i-XrB, (16) 
/mc^ \ c J 

in place of Eq. [71 Using the naive estimate Eq. [SI for B leads to a Lorentz factor of accelerated 
electrons, in place of Eq. [HI 



iGJe 



-3/4 /2 / GMrUp Y^ (I±\'^ 
\Gm) \h) \ ) \J^) 
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This may be rewritten in terms of fundamental constants, in analogy to Eq. [9l 

For a black hole of mass M ~ IO^Mq, as expected for AGN, electrons may be accelerated 
to energies ~ 10^° eV. For black holes of stellar mass the limiting electron energy is ~ 10^^ 
eV. If this energy is achieved in microquasars, they may produce (by Compton scattering) 
gamma-rays approaching this energy. 



4-2.2. Curvature radiation 

The extreme energies of Eqs. [161 [T7| and [18] are possible only in the limit in which Comp- 
ton drag is negligible, either because Lth is small or because is large (Eq. reflkn). If such 
extreme electron energies are achieved then even the curvature radiation may be powerful 
and energetic. Equating the accelerating power eEc ~ evB to the curvature radiation loss 
2e^C7^/(3r^) and using Eq. [6] for B yields a limiting Lorentz factor 

Gm) \h) [leJ [rj [m, 

. 4xio'f^y"V^)"Y^V'7^Y"\ (19) 

\1QGM j \10hJ \Le) \MqJ ' ^ ' 

where Vf. = e'^/niec'^ = 2.82 x 10^^'^ cm is the classical electron radius. Just as for pulsars, if 
Compton drag is negligible the electron energy will be limited by curvature radiation losses, 
and most of the power of electron acceleration will appear as curvature radiation. 

The characteristic frequency of curvature radiation is then 

1_ f lOGM \ 1/4 / /10GMq\ fmp\ ^'^ 

^ 27r V ) llO/iJ \Le) IOGMq V r,c^ ) \mj 



This frequency, typically hard X-rays for stellar-mass objects and very soft X-rays for AGN, 
is not a prominent feature of their spectra. This implies that curvature radiation is not 
usually the dominant electron energy loss mechanism in these objects, and supports the 
customary assumption that Compton loss is dominant. 
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4-2.3. Pair production 

Very energetic electrons lose more energy on a background radiation field by pair pro- 
duction (76" e~e~e'^) than by Compton scattering ((l3ll)). These authors find that this 
is the case only far into the Klein-Nishina regime, when 7e > bbOmeC^/huth- Using their 
results, we obtain a limiting Lorentz factor 



ik) ir m m & (^) ■ 

Of course, if this 'ye,pair exceeds the Goldreich- Julian limit Eq. [TSl then radiation drag is 
neghgible and the actual limit is 7e < 7Gje- 

Because of the sensitivity of 'je.pair to uth even a small amoung of soft radiation has a 
large effect on the limiting electron and positron energy, and hence on the photon energy of 
Compton gamma rays and of curvature radiation produced by these light leptons. This is 
the same difficulty encountered when Compton drag in the Klein-Nishina limit is dominant. 
Without a detailed geometrical and spectral model it is not possible to estimate 7e,pair even 
to order of magnitude. 



4.3. Proton acceleration 

If protons are accelerated their energy may be radiated as high energy gamma-rays 
following collisional pion production, either directly from 7r° decay or by Compton scattering 
of produced by decay of vr^. In the latter case neutrinos are also produced, with power and 
energy comparable to that of the gamma-rays. Alternatively, the proton kinetic energy may 
be degraded and coupled to electrons in a colhsionless shock, as is generally assumed to occur 
in gamma-ray bursts, and then radiated by Compton scattering or by synchrotron radiation. 
The proton kinetic energy may much exceed the values given by Eqs. U\ M and [9l and may 
approach the limiting energies given in §4.2.1. Gamma-rays (and neutrinos) following from 
pion production will have energies approaching a tenth of the proton energies. On the 
other hand, proton synchrotron radiation will occur at much lower energies than electron 
synchrotron radiation, and is generally negligible. 

If there were no pair production, the accelerated plasma would consist of protons (and 
nuclei) and electrons. Even in the presence of pair production, protons may be accelerated 
along with the positrons. This is important in intense sources of thermal radiation, such as 
BHXRB and AGN, because proton-photon scattering is neghgible below the pion production 



threshold. Even above threshold, the effective energy loss cross-section ((jl4j)) is a fraction 
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/ ~ 10 ^ of the Thomson cross-section. This permits accretion discs around black holes 
in AGN and BHXRB to be efficient proton accelerators (as has previously been discussed 



by (1301 : l26l : Il9l ) in other models). Very high energy gamma-rays may then result from 
photoproduction of 7r° by thermal radiation. Collisions of protons with nucleons make pions, 
leading to high energy radiation directly from vr'' decay and indirectly by Compton scattering 
of the very energetic produced by vr^ — >• ((19; 0))- Similarly, very high energy 



neutrinos are produced in the pion and muon decays. 

The energy loss length of a proton is, in analogy to Eq [5l 

1 Le rc^ 



i„ ~ -^^r. (22) 

Unlike Thomson scattering, this process has an energy threshold determined by the vr'^ rest 
mass of 135 MeV. If the thermal spectrum consists of visible light, the acceleration of protons 
is not restrained until 7p ~ 10^. The result analogous to Eq. [8]for the limiting Lorentz factor 
of protons slowed by photopion production in a thermal radiation field is 



^ \gm) [t-J (J \-^) ^ 



- ^ [wj [Togm) [fj (m) U) • ^''^ 

For an AGN with M ~ lO^M© this yields a limiting 7p ~ 10^, approximately the threshold 
at which 7r° photoproduction begins. Just as for pair production, the product of 7p and 
the black-body thermal photon energy is independent of M, so that Eq. [20] may be (barely) 
applicable at all M. 



4-3.1. Proton Radiation 

These possible ~ 10^'' eV protons in AGN could produce ~ 10^^ eV photons from pho- 
toproduced vr", but these are not observable at the greatest distances because of pair produc- 
tion on the microwave background radiation. However, at closer distances these gamma-rays 
may be observable. In BHXRB, perhaps including microquasars, the corresponding proton 
energies are ~ 10^^ eV, which may produce photons of ~ 10^^ eV. 
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4-3.2. Underluminous black holes as sources of UHE cosmic rays? 

If most of the thermal radiation is ineffective at slowing the protons because it is below 
the photopion threshold (in the protons' frame) or if Lth is small, then, just as for electrons in 
the Klein-Nishina case, it may be possible for protons to approach the limiting Lorentz factor 
Eq. [ini The low luminosity black hole powering Sgr A* in our Galactic Center demonstrates 
that black holes with Lth /Le ^ I may exist even in the presence of sufficient mass to support 
high accretion rates, and they would generally be undetectable at extragalactic distances. 
Ratios Lth/LE even lower than that of Sgr A* may occur. For parameters appropriate to 
supermassive black holes the proton energies of ~ 10^° eV, and energies of heavier nuclei 
greater by a factor of Z, may be sufficient to explain the highest energy cosmic rays. 



4.3.3. Stellar-mass black holes as TeV sources? 

Just as underluminous supermassive black holes may accelerate protons and nuclei to 
ultra-high energies, underluminous stellar-mass black holes may accelerate them to energies 
sufficient to be sources of TeV radiation (if the hadrons collide with other hadrons to make vr" 
or 71"^ whose decay electrons undergo Compton scattering). Such objects, resembling blazars 
(rather than isolated neutron stars, the more conventional explanation) might be the origin 
of the unidentified Galactic TeV sources ((1421)). 



5. Are there shocks in GRB? 

In most models of relativistic astronomical flows it is assumed that the fluids behave as 
if they were collisional, even though the single-particle collision lengths are generally orders 
of magnitude greater than the spatial extent of the fluids involved. Collective coUisionless 
processes (beam-driven electromagnetic or longitudinal two-stream instabilities) are assumed 
to lead to momentum transfer on much shorter scales, so that the flows can be described 
by the hydrodynamic equations. Kinematic constraints require some such interaction in 
order to convert the kinetic energy of relativistic motion to radiation. In particular, in a 
hydrodynamic supersonic flow the interaction takes the form of shocks. Such shocks are 
assumed in nearly all models of GRB. 

GRB are the only identifled objects near whose central engines X ^ 1. The observed 
radiation is produced at distances several orders of magnitude greater (typical estimates are 
010^^ cm), at which T <^ 1. Under these conditions nonthermal particle acceleration is 
possible, and in fact is required by their observed nonthermal spectra. Further development 



- 18 - 



of models requires investigation of the conditions under which thermal flows lead to the 
acceleration of nonthermal particles. Shock processes have been extensively investigated, 
but we should also ask if shocks (hydrodynamic discontinuities between two fluids each with 
thermodynamic equilibrium distribution particle functions, though generally transparent and 
with photon energy densities many orders of magnitude less than Planckian at the matter 
temperature) form at all. The kinematic constraints cannot be avoided, but it may be that 
the hydrodynamic assumption is invalid. It has been demonstrated ((0; S)) that the complex 
multipeaked temporal structure of GRB can only be explained, given the hydrodynamic 
assumption, by internal shocks. External shocks may be associated with, or even defined as. 



the origin of afterglows ( (|24|) ) . Yet elementary kinematics ( (l23l : l27l : |6|)) shows that, unless 
the ratios of Lorentz factors are very large, the energetic efficiency of internal shocks is low, 
generally no more than 10-20%. More efficient radiation would require the emission of a 
shell of low Lorentz factor followed by one of much higher Lorentz factor. It is implausible 
that this occurs regularly, but any deviation from this sequence (for example, a random 
distribution of Lorentz factors and of proper masses) would imply low efficiency. This would 
increase the energetic requirements of GRB and (perhaps worse) would raise the question 
of why the remaining energy does not appear as afterglows many times more energetic than 
GRB themselves. 



One possible escape from this "Efficiency Crisis" ((Il5l)) is to reject the assumption of 
hydrodynamic flow. Recall the reason why an external shock (on external matter initially at 
rest or moving non-relativistically) cannot explain complex substructure: After the external 
matter is hit by a shell of relativistic debris, it must be accelerated in a shock to relativistic 
speeds (Lorentz factor > 100) to explain the observed properties of GRB emission (spectrum 
and avoidance of gamma-gamma pair production). Then the external matter is no longer 
available as a static or non-relativistically moving target for the impact of a second shell, 
and the duration of emission of the first subpulse (from impact of the first shell) will overlap 
that of a second subpulse, even if a second shell is emitted cleanly separated from the first 
shell ((0;H)). 

This conclusion may be avoided if the first relativistic shell interacts only with a fraction 
of the matter at each point surrounding it, thus satisfying the kinematic constraints while 
violating the hydrodynamic assumption that the target behaves as a single fiuid with a 
single (relativistic) equilibrium particle distribution function. It is possible that the collective 
plasma processes that mediate the interaction couple, at any time, only to a small fraction 
of the matter (for example, to a narrow subrange of the velocity distribution) in the external 
medium. This is, at least qualitatively, consistent with the intermittent and fluctuating 
nature of most observed plasma-physical instabilities, both in the laboratory and in Nature. 
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This hypothesis suggests that as the density of target matter is gradually eroded by the 
passage through it of successive shells of relativistic debris the parameters of the accelerated 
particles and the properties of their radiation may also change. Although the details of these 
changes arc unknown (because the collective interactions between debris and target are not 
understood from first principles), this does suggest a progressive evolution of the charac- 
teristics of GRB emission through the pulse, consistent with the usual (but not universal) 
observations that the emission gradually softens. 

Unfortunately, this is not an unambiguous test of the hypothesis because the parameters 
and properties of the central engine arc also changing as its mass is gradually accreted on to 
a central black hole, and this can be an alternative explanation of any progressive changes 
through GRB pulses. The statistics of sub-pulses may be different in the internal shock and 
non-shock models, but in both cases depend on the unknown statistics of the activity of the 
central engine. 

6. Discussion 

It is apparent from comparing the results of §4.1 to those of §4.2 that the actual particle 
energies achieved depend sensitively on the flux and spectral distribution of any soft radiation 
in the acceleration region. The power and spectral distribution of radiation produced by 
interaction (Compton scattering, pion photoproduction) with this soft radiation also depend 
on the properties of the soft radiation, partly directly, and partly because it affects the 
distribution of energy of the energetic particles. A small amount of dense nonrelativistic 
plasma, producing a black body or thin bremsstrahlung fiux may have a large effect on 
the more energetic radiation of an accreting black hole, even though the thermal radiation 
is comparatively insignificant in power. The kinematic reason for this is similar to that 
encountered in the study of GRB external shocks, in which a proper mass E/V'^c^ at rest is 
sufficient to dissipate inelastically a kinetic energy E in matter moving with Lorentz factor 

r. 

In the present problem F may be as large as ~ 10^^ for at least a few of the energetic 
particles, and generally F ^ 1 by orders of magnitude, so that a tiny amount of isotropic 
thermal radiation (substituting for mass at rest) may have a large effect on the energetic 
particles if the interaction cross-sections are sufficient. This extraordinary sensitivity implies 
that predictive quantitative models are difficult to build. For example, for supermassive black 
holes possible objects range from classical quasars (with Lth ^ Le ^ lO*^*' crg/s) to sources 
of the most energetic cosmic rays with almost no electromagnetic luminosity at all. 
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Emergent spectra depend on such poorly understood and essentially unpredictable vari- 
ables as the flux and angular distribution of thermal radiation (and hence on the location 
and parameters of thermal gas) and on temporally fluctuating collective interactions between 
essentially collisionless interpenetrating streams of matter with relativistic relative velocities. 
Qualitatively, this sensitivity is consistent with the complex time structure observed in AGN, 
BHXRB and GRB, but makes it difficult to predict the emergent spectra. However, linear 
polarization is a general property of most nonthermal radiation processes (most familiarly, 
synchrotron radiation), and is predicted if the X-ray emission of BHXRB is nonthermal in 
some of their states. 

AGN may show QPO, just as suggested for GRB ((i23()). but with typical periods of 



order hours to days, depending on the black hole masses and angular momenta ((|40|)). 
No such QPO have been found in the extensive body of visible light data on AGN. This 
may be explained if the visible light is produced far from the central object or by thermal 
radiation from a nearly axisymmetric gas disc. It may be more fruitful (though more difficult 
observationally) to search for QPO in the energetic gamma-rays produced as particles are 
accelerated along magnetic field lines closer to the rotating disc. It is the magnetic field that 
would be expected to show the greatest deviation from axisymmetry, as in pulsars. 

The considerations of this paper have led to the prediction of two novel kinds of objects: 

1. Young, high field, rapidly rotating pulsars (with fi'^uj^ ^ 127Tc^Ichar ~ 10'''^ erg cm^/s^), 
that produce a black-body equilibrium pair gas wind rather than the few but very energetic 
particles produced by radio pulsars. Pulsars satisfying this condition on fi'^u^ may resemble 
a SGR in permanent outburst for a lifetime ~ 2u!^/{8nc^Ichar) ~ 3 x 10'"^^^ s^, where X is 
the moment of inertia. The minimum luminosity of such an object is Anr'^Ichar ~ 10^^ erg/s, 
so they would be detectable (following supernovae or quiet stellar collapse) as rapid (~ 1 
KHz) periodic X-ray sources at distances approaching that of the Virgo cluster (Ri 20 Mpc). 

2. Supermassive black holes with most of their luminosity in ultra high energy particles, and 
perhaps not recognizable as AGN. If these are the source of UHE cosmic rays, their arrival 
directions will be clumped (subject to propagation effects) and perhaps correlated with point 
sources of UHE neutrinos and gamma-rays. 



I thank T. Piran, M. A. Ruderman and R. Sari for discussions, Washington University 
for the grant of sabbatical leave, and the Hebrew University for hospitality and a Forchheimer 
Fellowship. 



REFERENCES 

Barraud, C, et al. 2003, A&A 400, 1021 
Blandford, R. D. 1976, MNRAS 176, 465 

Blumenthal, G. R. & Gould, R. J. 1970, Rev. Mod. Phys. 42, 237 

Dar, A. & Laor, A. 1997, ApJ 478, 5 

Eggum, G. E., Coroniti, F. V. & Katz, J. I. 1988, ApJ 330, 142 



Fan, Y. & Piran, T. 2006, MNRAS 369, 197 ( astro-ph/0601054D 
Fenimore, E. E., Madras, C. D. & Nayakshin, S. 1996, ApJ 473, 998 
Frail, D. A., Waxman, E. & Kulkarni, S. R. 2000, ApJ 537, 191 
Frail, D. A., et al. 2001, ApJ 562, L55 
Gaidos, J. A., et al. 1996, Nature 383, 319 
Galama, T. J., et al. 1998, ApJ 500, L97 

Gnedin, Yu. N., Borisov, N. V., Natsvlishvili, T. M., Piotrovich, M. Yu. & Silant'ev, N. A. 
2003, |astro-pli/0304 158l 

Goldreich, P. & Julian, W. H. 1969, ApJ 157, 869 

Greisen, K. 1966, Phys. Rev. Lett. 16, 748 



loka, K., Toma, K., Yamazaki, R. & Nakamura, T. 2006, A&A 458, 7 (astro-ph/0511749D 

Jones, F. C. 1965, Phys. Rev. 137B, 407 

Katz, J. I. 1977, ApJ 215, 265 

Katz, J. I. 1982, ApJ 260, 371 

Katz, J. I. 1991, ApJ 367, 407 

Katz, J. I. 1994a, ApJ 437, 724 

Katz, J. I. 1994b, ApJ 432, L107 

Katz, J. I. 1996, ApJ 463, 305 



- 22 - 

Katz, J. I. 1997, ApJ 490, 633 (| astro-ph/970n76l ) 
Katz, J. I. 1998, A&AS 138, 531 
Katz, J. I. & Salpeter, E. E. 1974, ApJ 193, 429 
Kazanas, D. & Ellison, D. C. 1986, ApJ 304, 178 
Kobayashi, S. & Piran, T. 1997, ApJ 490, 92 

Krolik, J. H. 1999, Active Galactic Nuclei (Princeton: Princeton U. Press) 
Lovelace, R. V. E. 1976, Nature 262, 649 
Lynden-Bell, D. 1969, Nature 223, 690 

Mastichiadis, A., Marscher, A. P. & Brecher, K. 1986, ApJ 300, 178 

Nemiroff, R. J., et al. 1994, ApJ 423, 432 

Piran, T. & Shaham, J. 1975, Nature 256, 112 

Piran, T. & Shaham, J. 1977, ApJ 214, 268 

Piran, T., Shemi, A. & Narayan, R. 1993, MNRAS 263, 861 

Press, W. H. & Rybicki, G. B. 1997, Astronomical Time Series, eds. E. Leibowitz, D. Maoz 
& A. Sternberg (Dordrecht: Kluwer) p. 61 

Quinn, J., et al. 1996, ApJ 456, L83 

Ruderman, M. A. & Sutherland, P. G. 1975, ApJ 196, 51 

Salpeter, E. E. 1964, ApJ 140, 796 

Sari, R. & Piran, T. 1997, ApJ 485, 270 

Shapiro, S. L. & Teukolsky, S. A. 1983, Black Holes, White Dwarfs and Neutron Stars (New 
York: Wiley), 361 

Sturrock, P. A. 1971, ApJ 164, 529 

Ubertini, P. 2005, Experimental Astronomy 20, 75 (|astro-ph/0605228) 

Weiskopf, M. C., Sutherland, P. G., Katz, J. I. & Canizares, C. R. 1978, ApJ 223, L17 

Zdziarski, A. A., Poutanen, J., Paciesas, W. S. & Wen, L. 2002, ApJ 578, 357 



-23- 



This preprint was prepared with the A AS macros v5.2. 



